A correction method for the pair density to get close to the ground state one 
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We present a correction method for the pair density (PD) to get close to the ground state one. 
The PD is corrected to be a variationally-best PD within the search region that is extended by 
adding the uniformly-scaled PDs to its elements. The corrected PD is kept A'^-representable and 
satisfies the virial relation rigorously. The validity of the present method is confirmed by numerical 
calculations of neon atom. It is shown that the root-mean-square error of the electron-electron 
interaction and external potential energies, which is a good benchmark for the error of the PD, is 
reduced by 69.7% without additional heavy calculations. 
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The electron correlation is one of main topics in the 
fields of the atomic, molecular and condensed mat- 
ter physics. The physical quantity that directly ex- 
presses the electron correlation is the diagonal element 
of the second-order reduced density matrix, which is so- 
called the pair density (PD). The PD gives not only 
the electron density but also expectation values of ar- 
bitrary two-particle operators flLISI. Therefore, the PD 
functional theory (PDFT) j4i-i33i] is one of the promis- 
ing first-principle methods for describing the electron 
correlation beyond the conventional density functional 
theorvfs^. [ssj . In this letter, we present a correction 
method for the PD to get close to the ground state one, in 
which the search region of PDs is substantially extended. 
This method can become a powerful tool to overcome the 
well-known problems of the PDFT. 

There exist two kinds of well-known problems in de- 
veloping the PDFT. One is to have to develop the ap- 
proximate form of the kinetic energy functional since the 
kinetic energy cannot be expressed by the PD alone [3- 
Some approximation is needed for the kinetic 
energy functional of the PDFTjlB-il. The other prob- 
lem is related to the iV-representability problem of the 
Although we have to search the ground 
state PD within the set of the iV-representable PDs (2nd 
Hohenberg-Kohn (HK) theorem of the PDFT) , the nec- 
essary and sufScient conditions for the iV-representability 
of the PD have not yet been known in a practical form[l|- 
1 [23423. 

In this letter, concerning the first problem we propose 
an approximate functional of the kinetic energy along 
the strategy shown in the previous paper [13. Namely, 
an approximate functional is developed by imposing two 
kinds of sum rules for the kinetic energy functional as the 
restrictive conditions. Although the previous approxi- 
mate functionals have been developed by imposing exces- 
sively the restrictive conditions on them[19], the present 
one is devised with removing such excessive restrictive 
conditions 13, 33| • We have 
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7(^)(rr';rr')drdr', 
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where j^-^^rr'; rr') denotes the PD, and where K and K' 
are arbitrary constants that are determined later. Equa- 
tion ll]) satisfies the scaling property of the kinetic energy 
functional and is consistent with the HK theorem of the 
PDFT[ll|3i|. The former fact will be significant later 
in discussing a correction method. 

Concerning the second problem, we have recently 
attempted to extend the variational search region of 
PDs 17, l3| beyond that of the effective initial theory 13, 
ITsj . However, it has been eagerly anticipated that such 
a search region is further extended with keeping the A'^- 
representability of the PD. In this letter, we shall extend 
the search region by adding another set of PDs to the 
elements of the search region. This is precisely our cor- 
rection method, and is a main subject of this letter. Let 
us show this method in the following three steps (i), (ii) 
and (iii). 

(i) First of all, we shall show that the kinetic energy 
functional which is consistent with the scaling property 
satisfies the virial relation exactly. Here suppose that 
the search region of PDs fully covers a set of the A^- 
representable PDs. Applying the scalings of the electron 
coordinates, 7*^^'' (rr'; rr') is transformed into 



7f^(rr';rr') = AS^^H^rAr'; ArAr'), 
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where A is the scaling parameter. 73^ (rr'; rr') is the PD 
scaled from 7(^'(rr'; rr'), which is called the scaled PD. 
The kinetic energy functional that is consistent with the 
scaling property satisfies the following relationjsl [Til [l3| : 



rr (2) 



v(2) 



(3) 



2 



Similarly, the electron-electron interaction and external 
potential energies to the scaled PD are, respectively, 
given by 
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-j7(2)(rr';rr')drdr', (5) 



where t'est(r) denotes the external potential and where 
N is the number of electrons. Due to the 2nd HK 
theorem of the PDFT, the total energy functional 
E [7(2)] T [7(2)] + W [7(2)] + V [7(2)]) is minimum 
(2) 

at 7q that corresponds to the variationally-best PD 
among all of the iV-representable PDs. Using the scaled 
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PD, the theorem can be written as dE 

(2) (2) 

0. Here 7q 3^ is the PD scaled from 7g , and is defined 
similarly to Eq. ©. Substituting Eqs. ©, ^ and © 
into this theorem, we have 
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Equation (jS]) is the virial relation of the PDFT. It should 

be noted that even though the kinetic energy functional 

that satisfies Eq. ([3]) is not an exact but an approximate 

form, Eq. (j6]) exactly holds if the search region fully 

covers the set of the A^-representable PDs. 

(ii) Next, we shall consider the case where the search 

region not fully but partially covers a set of the iV- 

representable PDs. Suppose that such a partial search 

region is denoted as Ci (See, Fig. 1), and that the total 

~ (2) 

energy functional is minimum at 7q that corresponds to 

the variationally-best PD within the set Ci. Due to the 

incomplete cover of iV-representable PDs by the set Ci , it 

(2) 

could be that not all of the scaled PDs 7q 3^ are included 
in the set Ci. We denote the set of 7q ^ as C2 (Fig. 1). 



In this case, the total energy functional E 
not always take the minimum at A = 1 



7o'i] does 
Namely, we 



have dE 



^(2) 
'0, A 



'd\ 



7^ 0. Like the derivation of Eq. 



([6|), this inequality leads to the relation such that the 
left-hand side of Eq. ([6|) is not equal to the right-hand 
side. Thus, if the search region does not fully cover the 
set of the iV-representable PDs, the virial relation does 
not hold even though the kinetic energy functional sat- 
isfies Eq. ([3]). This fact can be used as a criterion of 
whether the search region fully covers the set of the A^- 
representable PDs or not. Here let us define the virial 
ratio which indicates to what extent the virial relation 
holds: 
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Using this ratio, the above-mentioned criterion can be 
restated as follows; the deviation of from the value 
-2.0 means the insufficiency of the search region of PDs 
if the kinetic energy functional satisfies Eq. /f^)[37|. 

(iii) When Ry deviates from -2.0, a correction method 
that is related to the extension of the search region is 
desired. The key point of the present correction method 
is the extension of the search region by adding the scaled 
PDs to the search region. Figure 1 truly shows the re- 
lation between the original search region Ci and newly- 
added search region C2. Both sets Ci and C2 are sub- 
sets of the set of the iV-representable PDs. The above- 



mentioned inequality dE 
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^ means that 



the set C2 could possibly include the PD that takes a 
total energy lower than what the best PD within the set 

(2) 

( ' . i.e., 7g , takes. Such a PD exists in C2 n Ci and is 

(2) ~ 

denoted as 7q j^, where Ci is the complementary set of Ci 



(Fig. 1). Then, the inequality dE 
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0. Substituting Eqs. 



([3|), Q and ^ into this equation, we obtain 
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If we consider the isolated atomic system, Eq. ([5]) is 
easily rewritten as 
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From Eq. ([5]) or we can determine the value of 

(2) 

A. The scaled PD 7^ ^ with A thus determined is the 

variationally-best one within the union of two sets, i.e., 

Ci U C2 . That is to say, a correction of the ground state 
~ (2) ~ (2) 

PD from 7q to 7q ^ is accomplished by extending the 
search region from Ci to Ci U C2. It is obvious that 

(2) 

the corrected PD 7^ remains A^-representable. Fur- 
thermore, it is easily confirmed that the virial relation 

(2) 

exactly holds for the corrected PD 7^ j^. 

The procedure of this correction method is as follows. 
First we prepare the approximate form of the kinetic en- 
ergy functional that satisfies Eq. ([3]), and determine the 
search region of PDs. Using this kinetic energy func- 
tional, we calculate 7q . Substituting 7g into Eq. ^ 

or ([9]), the scaling parameter A is obtained. This A im- 

(2) 

mediately leads to the corrected PD 7^ 3^ by using Eq. 
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FIG. 1: The relations between various sets of PDs. The set 
Ci covers a set of the A'^-representable PDs not fully but par- 
tially. The set C2 consists of the uniformly-scaled PDs. The 
variationally-best PD within Ci, and that within Ci UC2, are 
denoted by 7q and 7^ ^ , respectively. 



([2]). Of course, the kinetic energy, electron-electron inter- 

~ (2) 

action energy and external potential energy for 7q ^ can 
be easily calculated from Eqs. ([3|), (HI and ([5]), respec- 
tively. One of striking points of this correction method 
is that the heavy calculation task is not needed for the 
corrections of the PD and related quantities. 

In order to check the validity of the present method, 
we perform numerical calculations for neutral neon atom. 
We utilize the computational PD functional scheme [l9j| 
so as to obtain t-q^'. According to this scheme fioj. we 
adopt as the search region the set of PDs that are calcu- 
lated from the correlated wave functions. As the corre- 
lated wave functions, we take the linear combination of 
the ground state and doubly-excited Slater determinants 
(SDs) that consist of eigenfunctions of the effective ini- 



tial scheme 3, 



18|. Among all doubly-excited SDs, we 
choose the SDs that have non-negligible contributions to 
the resultant PD. A total of 2,861 SDs and 12,275 SDs, 
the constituent eigenfunctions of which have the princi- 
pal quantum number up to 6 and 11, respectively, are 
chosen as basis functions in constructing PDs. 

As the set Ci, we take the set of PDs that are con- 
structed by 2,861 SDs. First, by means of the computa- 
tional PD functional scheme[l^, we calculate 7q^^ that is 
the variationally-best PD within the set Ci. Then, the 
corrected PD 7^^]^ is calculated by using Eqs. ^ and (0). 
We also prepare the search region that is extended with 
increasing the number of basis SDs (NSD) from 2,861 to 
12,275. Hereafter, we denote this search region of PDs as 
C'l . For comparison, the variationally-best PD within the 
set C'l is also calculated by means of the computational 
PD functional schemellOl. 



Using 7^- 
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and 7q^ 



A' 



we calculate the errors of the 



electron-electron interaction energy, external potential 
energy, kinetic energy and virial ratio that are denoted 
by AVF, AF, AT and Ai?^, respectively. The root- 
mean-square error (RMSE) of the electron-electron inter- 
action and external potential energies is also calculated 
by ^{(AW^)2 + (AV^)2}/2. The RMSE is considered as 
a good benchmark for to what extent the resultant PD 
is close to the correct ground state PD. This is because 
accuracies of both the electron-electron interaction and 
external potential energies are dependent only on that of 
the PD. 

The above-mentioned calculations are performed with 
changing the values of K and K' that appear in Eq. ([T|) . 
These values are determined by requiring the approxi- 
mate functional to have two desirable features. One is 
that the approximate functional reproduces the Hartree- 
Fock kinetic energy if the PD coincides with the Hartree- 
Fock PD. The other is that the RMSE of the electron- 
electron interaction and external potential energies is 
minimized with respect to parameters. This determi- 
nation process is implemented in individual cases before 
and after the correction. 

Calculation results are shown in Fig. 2. From this 
figure, we find the following points. 

(I) It is found in Fig. 2 that the present method re- 
duces AVF and AF by 89.7% and 66.3%, respectively 
(reductions from case a to c for "AVF" and "Ay, re- 
spectively, in Fig. 2). On the other hand, only 1.3% and 
1.1% reductions are made via increasing the NSD from 
2,861 to 12,275 in AM^ and AF, respectively (reduction 
from case a to 6 for "AVF" and "AF", respectively). 
Also, as shown in Fig. 2, the RMSE is remarkably re- 
duced by 69.7% (reduction from case a to c) , while it 
is improved only 1.2% by increasing the NSD (reduction 
from case a to 6). Thus, the reduction rate of the present 
method is much larger than that of the extension method 
with increasing the NSD. This leads to that the search 
region of PDs is more effectively extended by adding the 
uniformly-scaled PDs than by increasing the NSD. Fur- 
thermore, it should be noticed that the present method 
is feasible without additional heavy calculations. 

(II) It is found in Fig. 2 that AT is also reduced much 
more effectively by the present method than by increas- 
ing the NSD (80.9% and 1.2% reductions for the former 
and latter, respectively). In general, the accuracy of the 
kinetic energy is dependent on both the appropriateness 
of the approximate form of the functional and that of 
the resultant PD. Judging from the results of the RMSE, 
/S.W and ISV (Fig. 2), the resultant PD seems to be im- 
proved to be close to the ground state one. Therefore, 
we can deduce that the approximate form itself, which is 
given by Eq. (1), would be also sound. 

(III) It is found in Fig. 2 that the virial ratios Ry 
before the correction deviate from the correct value -2.0. 
Using the criterion on the search region, these deviations 
mean that not only the search region Ci but also the 
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FIG. 2: The calculation results of the AW, AV, RMSE, AT 
and ARv In this figure, the search region Ci, C'l and Ci UC2 
are simply denoted by the symbols " a" , " b" and " c" , respec- 
tively. 



search region C( are insufficient since the approximate 
functional Eq. ([IJ is consistent with Eq. That 
is to say, even though more than 12,000 SDs are used 
in constructing PDs, the search region is not extended 
effectively. On the other hand, Ry after the correction is 
rigorously equal to -2.0, as it should be. These tendencies 
seem to be consistent with those mentioned in the above 
(I) and (II). 

Thus, the present method improves not only AW, AV 
and RMSE but also AT and ARy quite substantially and 
effectively. We can therefore say that the PD is corrected 
to be close to the ground-state one appropriately via the 
present method. 

In conclusion, we summarize the features of the correc- 
tion method proposed here. The most distinctive feature 
is that 

(1) the search region is extended by adding a set of the 
scaled PDs to elements of the search region. 

The resultant PD that is corrected by this method pos- 
sesses the following features: 

(2) the corrected PD is A^-representable, 

(3) the corrected PD satisfies the virial relation exactly. 
The validity of the present method is successfully con- 
firmed by numerical calculations of neon atomic system. 
It is shown that 

(4) not only the RMSE but also the errors of the kinetic 
energy, electron-electron interaction energy and external 
potential energy are individually all reduced definitely, 

(5) the correction is enough effective even in the small 
size calculations. 

This correction method provides an ingenious way to 
extend the search region of PDs efficiently. Due to the 
computational easiness of the correction, this method 
will be extensively applicable to larger systems such as 
molecules, clusters and solids. 
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